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ABSTRACT 

LeBlanc, Philip Norman. M.S. Purdue University, August , 1972. Land Use 
Classification Utilizing Remote Multispeetral Scanner Data and Computer 
Analysis Techniques. Major Professors: Christian J. Johannsen and 

Joseph E. Yahner. 

This research was designed to study the ability of present automatic 
computer analysis techniques with the use of multispeetral scanner data to 
differentiate land use categories represented in a complex urban scene and 
in a selected flightline. An airborne multispeetral scanner was used to 
collect the visible and reflective infrared data. 

A small subdivision near Lafayette, Indiana was selected as the 
test site for the urban land use study. Multispeetral scanner data were 
collected over the subdivision on May 1, 1970 from an altitude of 915 
meters. The data were collected in twelve wavelength bands from 0.1*0 
to 1.00 micrometers by the scanner. 

The results indicated that computer analysis of multispeetral data 
can be very accurate in classifying and estimating the natural and man-made 
materials that characterize land uses in an urban scene. 

A 1.6 km. wide and 16 km. long flightline located in Sullivan County, 
Indiana, which represented most major land use categories, was selected for 
analysis. Multispeetral scanner data were collected on three flights from 
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an altitude of 1,500 meters. Energy in twelve wavelength bands from 
0.46 to 11.70 micrometers was recored by the scanner. 

A new, more objective approach to computer training w as developed 
for analysis of the three dates of data. Emphasis was placed on the 
standardization of a procedure for analysis of data. The procedure 
offered faster and consistantly good duplication of attained results. 

The results' indicated an ability for automatic computer analysis 
of remotely sensed multispectral scanner data to characterize and map 
land use categories within the test area. Additionally, results in- 
dicated an alteration of the data analysis procedure and land use class- 


ification scheme 
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INTRODUCTION 


The need for current information on major land use, acreage, and 
distribution of crops is essential for efficient management of resources 
in both the developed and developing countries of the world. Histori- 
cally, the collection, compilation and interpretation of data on use of 
land in the United States has proved invaluable in the study of present 
land-resource problems. 

In the past, land use data have been compiled by census interview 
or by field mapping. In some countries, major classes of land use are 
mapped in the field and published, usually on a small scale. For most 
countries, land use inventory data consist of data compiled by census 
from personal interview, mail questionnaires, study of sample areas, or 
some combination of these means. Generally, these methods take a re- 
latively long time because the number of trained scientists and workers 
are limited. These factors account for the relatively long intervals 
between census projects in many countries. 

The increasing need for land use inventory in many developed and 
some of the developing countries has already pressed the use of one 
method of remote sensing. This familiar technique 13 the use of black 
and white aerial photography. Air photos have been used to make mans 
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showing the distribution of major land uses and specific crops in many 
countries. Measurements from these airphoto maps provide data on 
acreage of categories of land use. Comparison between recent and older 
photos can provide information regarding change such as expansion of 
urban development , revision of farmland to forest , or development of 
new farmland by clearing, drainage, or irrigation. 

Currently, the competition for use of land is becoming more 
intense. Urban development, need for more recreation areas (particu- 
larly near large urban areas) and the preservation of wildlife habitats 
are matters of great interest to those concerned with land resources. 

In the future, it is clear that current, accurate and concise informa- 
tion will be of great importance and need in the careful examination 
and utilization of our resources. 

Newer, more sophisticated remote sensing systems have the poten- 
tial to offer significant improvement in rural and urban land use eva- 
luation. Computer analysis of multispectral scanner data is one such 
new am in the field of remote sensing. This project is concerned with 
determining the ability of present automatic computer techniques 
employing multispectral scanner data to obtain useable land use infor- 
mation. 

The first study is a preliminary investigation into the uses of 
remote sensing techniques for the identification of land use3 that 
characterize an urban scene. The second study involves the use of a 
new multispectral analysis technique for the identification of land use 
classes in an agricultural scene. The results from this project should 
further define the role of remote sensing as a tool in resource evalua- 


tion. 
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CHAPTER 1 
LITERATURE REVIEW 


The Role of Remote Sensing in Land Use 
Classi f i cation 

The concept of multispectral remote sensing is that "plants and 
materials located at the surface of the earth possess the fun drum :ntal 
property of reflecting or emitting electromagnetic energy differentially 
over the spectrum" (Tanguay, 1969a). In the visible portion of the 
electromagnetic spectrum these particular properties are seen in the 
form of color. These properties appear to be characteristic of each 
material and are due to the fundamental properties of energy absorption, 
emission and surface reflectance. These characteristics will be further 
discussed in following sections. 

Hoffer (1967a) has defined multispectral remote sensing as: 

"...the detection and recording, from a distance 
or remote location, of reflected or emitted 
electromagnetic radiation in many discrete 
relatively narrow spectral bands between 0.3 
to 14 micromotors wavelength..." 


Remote sensing information can be collected by a variety of instru- 
ments mounted near or above the surface of the earth in airplanes and 
satellites. With these instruments, energy reflected from a "target" 
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can be collected and recorded for processing and analysis. Holter 
(1970) and LARS (1968) have described in detail the instrumentation, 
data collection and data processing systems currently used in the field 
of remote sensing with multispectral scanner data. 

The merits of remote sensing as a potentially valuable tool in 
gathering and updating land U3e and resource information has been 
reported by Shelton (1967) and Johannsen and Baumgardner (i960). Aside 
from its data acquisition and handling characteristics, the value of 
remote sensing in a land use classification system depends on its 
ability to accurately differentiate materials such as rock, bare soil, 
crop species, trees, grass, and water. The basic assumption underlying 
this project is that remote sensing, utilizing multispectral scanner 
data and automatic computer analysis techniques, has achieved sufficient 
accuracy and reliability in the identification of land surface materials 
to warrant its application to land use classification. The literature 
reports successful application of remote sensing in the areas of vege- 
tation mapping, soils mapping and geologic studies. 

In soils mapping, studies by Kristof and Zachary (1971), Stoner and 
Horvath (1971V West (1972), have shown that spectral maps of surface 
soils can he produced using multispectral data and automatic computer 
analysis techniques which compared favorably with soils maps prepared 
by conventional soil survey techniques. Mathews et al. (1971), using 
multispectral data, were able to classify soil parent materials that 
compared well with the boundaries identifiable on aerial photography. 

It has been shown (Al-Abbas et al. , 1972) that gross textural differences 
in surfo.ee soils can be mapped using remote sensing techniques. 
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Tanguay ( 1969 b), in research studying the feasibility of engineering 
soils mapping, reported identifying and mapping areas that indicated 
"the distribution of main soils classes, drainage features, wet zones, 
muck pockets, and bare rock areas." In another study, Baumgardner et 
al. (1970) have shown that computer analysis techniques can differen- 
tiate five different ranges of organic matter content for mineral soil3. 

Recent geologic studies have reported classifying various geologic 
materials. Srnedes et al. (1969) have reported using multispectral data 
and automatic computer analysis techniques to map terrain in Yellowstone 
National Park. They showed that nine classes (bedrock exposures, talus, 
vegetated rock rubble, glacial kame, glacial till, forest, bog, surface 
water and shadows) could be mapped with good accuracy. Using multi- 
spectral scanner data collected between 6.7 and 13. U micrometers and 
computer analysis techniques; Lyon et al, (IO 69 ) were able to identify 
quartz beach sand at a Texas site and recognize and map granite and 
andesite rock at a California test site. In preliminary work with 
resource inventory and survey, Mallon (l97l)> using scanner data collected 
in the 8 to 1^ micrometer range of the electromagnetic spectrum, inves- 
tigated automatic computer mapping of manganese, bauxite and chromate ore 
materials. He reported that his studies show "an appearance of a family 
of spectral curves" associated with each material. 

In vegetation mapping studies, correct identification of 8l to 99 
percent were reported in mapping six cover-type classes (row crops, 
wheat, oats, pasture, hay and trees) using multispectral scanner data 
(LARS, 1970). Hoffer ( 1967 b) reported good results in classifying land 
surfaces into thi'ce basic categories of green vegetation, bare soil and 


water. 


The literature indicates that basic land cover mapping utilizing 
multispectral scanner data and presently developed computer analysis 
techniques is sufficiently advanced to provide accurate identification 
and mapping of land use characteristics. It is with this assumption 
that this study was undertaken. 

Energy Phenomena 

The electromagnetic spectrum has no definite upper or lower limit 
it can be stated, however, that it extends from radio frequencies , less 
than one Hertz (cycles per second) to cosmic rays greater than 10 26 
Hertz (Halliday and Resnick, 1962). Hoffer and Johannsen (1969) illus- 
trated a portion of the electromagnetic spectrum which is of importance 
to remote sensing (Figure 1) and reviewed various types of sensors used 
in the different wavelength regions. 

Optical-mechanical scanner systems effectively operate within the 
electromagnetic spectrum from approximately 0.3 to 15.0 micrometers. 
This portion of the spectrum can be divided essentially into two parts : 
the reflective and the emissive regions (approximately 0.3 to 3-0, and 
3.0 to 15.0 micrometers, respectively). The reflective region can be 
further subdivided into the ultraviolet. Visible and infrared regions 
(approximately 0.004 to 0.38, 0.38 to 0.68, and 0.68 to 3.0 micrometers 
respectively). 

A remote sensor receives energy emanating from a target as a 
function of the illumination and the reflective and emissive properties 
of the target. The illumination source for the earth is primarily the 
sun; although, illumination may be received from the moon, stars, or 
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intense artificial sources such as a laser or mercury- arc lamp. The 
earth's atmosphere and local weather changes can .adversely affect the 
spectral character of the sunlight (Gates, 19 6 5 ; and Johannscn , 19o9). 

The attenuation effects by the atmosphere on the sun's energy profile 
must be considered not only from the source to the target, but also the 
energy returning from the target to the remote sensor. Within the elec- 
tromagnetic spectrum there are regions or "windows" where the attenu- 
ating effects of the atmosphere are relatively hindered. Within these 
regions spectral channels are selected for energy measurements from a 
target to a remote sensor. 

Radiated energy is emitted by all objects above the temperature of 
absolute zero as a function of their temperature and emissivity 
(Gates, 1962). Basic data on emissivity are given by Myers and Allen 
(1968) and Kumar (1972). Various factors influence soils and plant 
temperatures. Gates (.1963), Waggoner and Shaw (.1951), Ansari and Loomis 
(1959) and Kryuchkow (1961) discussed the factors of solar radiation, 
air temperature, humidity, wind speed, available soils moisture, leaf 
structure and differences in plant species as the important factors which 
affect plant temperature. The primary factors relating to soil tempera- 
ture are soil particle size, soil moisture, soil air humidity and soil 
porosity (Myers, 1970). 

The reflectance function varies considerably from one object to 
another and for each object within :i specific category. Gates (1970) , 
Gausman et al. (.19^9) and McClure ( l'X'9 ) discussed many of the physical 
and physiological properties of plants which explain some of their re- 
flectance character! y. t i cs . The properties mentioned were leaf geometry, 
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morphology, physiology, chemistry, soil silo and climate. Coil reflec- 
tance characteristics are influenced by moisture content, particle size, 
structure, surface roughness and organic matter content (Bowers and 

Hanks, 1965). 

From this discussion it can be seen that objects reflect and emit 
energy differentially over the electromagnetic spectrum. It is the 
concept of multispectral analysis that in some channel or channels in the 
electromagnetic spectrum classes or categories of materials will be 
distinguishable and will be detectable with multispectral analysis tech- 
niques. It must be kept in mind that each factor mentioned that influ- 
ences the reflectivity, emissivity or illumination of an object can also 
contribute considerable variability to the response of a given material. 
It is for this reason that frequently there will exist several spectral 
subclasses to characterize one cl iss of material. 

Land Use Clas sification Schemes 

Because there are several concepts about what "land" is 
(Barlowe, 1958 ), there exists a similar range in ideas about the concept 
of "land use." In an effort to reach a general understanding of the 
definition of the term "land use" the Committee on Land Use Statistics 
grappled at length with this problem (Clawson and Stewart, 1965)- The 
committee concluded that the term "land use" should refer "to man's 
activities on land which are directly related to the land." It was 
further noted that phrases such as "human uses of land" or "human acti- 
vity on land" might help in the clarifying of the concept of land use; 
it is the activity dimension that constitutes the central focus in 
defining land use. 
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The purpose of this section is to review some of the presently 
accepted land use schemes in on effort to select or formulate a suitable 
scheme that would be compatible with remote multispectrally sensed 
scanner data. Four land use schemes from the Conference on Land Use 
Information and Classification (1971) were selected for review; a more 
detailed outline format of each scheme is included in the Appendix, 

Table B. 

Three of these classification schemes, the Standard Land Use 
Coding Manual, the Canadian Land Inventory, and a scheme developed by 
the Association of American Geographers, were developed for use on a 
country-wide or national basis, while the fourth scheme, the f/ew York 
State Land Use scheme, was devised for use only in the state of flew York. 

Tiie three national schemes differ significantly in the level of 
inventory detail contemplated in their use. The Standard Land Use Coding 
Manual ("A Standard System for Identifying and Coding Land Use Activities," 
1965) provided a four-digit categorization of land use which was developed 
mainly for identifying and coding (rather than mapping) urban land use 
C activities and adjacent situations in the United States. This classifi- 
cation scheme was not designed specifically for use with air photo-inter- 
pretation or other remote sensing technques. Ground observations and 
enumeration must provide much of the information necessary to classify 
land use with this scheme when used in urban areas. The scheme has nine 
first-order categories and ninety-nine second-order categories. 

The Canadian Land Inventory ("A Guide to the Classification of Land 
Use," 1970) employed essentially a two-digit or two level land use classi- 
fication scheme for use throughout the more densely populated portions of 
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Canada. Maps of present land use were complied at a scale of 1:50,000 
using this scheme, which has six major or first-order categories, twelve 
second-order categories, and two third-order categories. The scope and 
objectives of the Canadian inventory necessitated a much more generalized 
scheme for classifying land use than the situations for which the 
Standard Land Use Coding Manual was primarily prepared. 

The third national scheme was developed from a study by the 
Association of American Geographers (Anderson, 1971 a) for use with 
orbital imagery for making land use maps ranging in scale from 1:250,000 
to 1:2,500,000. This classification scheme vus developed and used in a 
pilot study in the Phoenix, Arizona area to test the capabilities of the 
devised scheme for use with conventional color and color infrared imagery 
taken from the Apollo 9 satellite and from high altitude aircraft (50,000 
to 60,000 feet). The imagery was supplemented by black and white photo 
mosaics at a scale of 1:62,000. The scheme has seven first-order cate- 
gories and twenty third-order categories. Fourth or even fifth-order 
categories can be added for use with larger scale imagery or for use with 
ground observations and enumeration. 

The fourth scheme developed specifically for the state of Hew York 
is quite adaptable for use elsewhere. The New York land inventory classi- 
fication scheme ("A Review of the New York State Land Use and Natural 
Resource Inventory," 1970) employed essentially a three-digit classifica- 
tion scheme which was developed to accomodate the recognition of areal and 
point data, i.e., cropland and farmsteads. The scheme has eleven first- 
order categories, forty-eight second-order categories and seventy-five 
third-order "point data" categories. Aerial photography at a scale 
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1 ,000 constituted the main land use information source for the hew 

York state inventory. 

A set of working criteria against which to evaluate land use clas- 
sification schemes for use with orbital and other high altitude imagery 
was presented by Anderson (1971b). It proved to be valuable in the 
final evaluation of the selected classification scheme used in this 
study. These criteria are: 

(1) A minimum level of accuracy of about 85 to 90 percent 
or better should be approached in the interpretation 
of the .imagery being used. 

(2) A well-balanced reliability of interpretation for the 
several categories included in the classification 
scheme should be attained. Closely related to the 
requirement of a minimum overall level of accuracy 

is the matter of varying levels of accuracy which can 
be attained for the several categories of the classi- 
ficati j .\ being used. 

(.3) Repeatable or repetitive results should be obtainable 
from one interpreter to another and from one time of 
sensing to another. 

( 4 ) The classification scheme should be useable or adapt- 
able for use over an extensive area. 

( 5 ) The categorization used .in the classification scheme 
should permit vegetative and other cover types to be 
used as surrogates for activity-oriented categories 
wherever possible. 

(6) The classification scheme should be suitable for use 
with imagery taken at different times during the year. 

( 7 ) The classification scheme should permit effective use 
of sub-categories that can he obtained from ground 
surveys or from the use of imagery available at a 
larger scale or with the use of color photography. 

(8) A need to collapse the categories of the classification 
scheme into a smaller number of categories mu3t be 
recognized. 
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(9) The classification scheme should recognize the multiple 
use aspects of land use whenever possible. 

The above outlined criteria could not be met by the classification 
scheme developed for use in this study. As Anderson (1971) has stated, 
even some of the presently used classification schemes do not meet 
satisfactorily the above criteria. In the future, it may be desirable 
to augment these criteria as the role of remote sensing in a land clas- 
sification system is more fully defined. 
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CHAPTER II 

TRAINING FIELD APPROACH TO CLASSIFYING AN 
URBAN SCENE — A PRELIMINARY STUDY 

Int roduction 

The uso of remote sensing techniques in the rapid gathering, pro- 
cessing and interpretation of data obtained over an urban area is of 
interest to many people. Quantitative information from a remote sensing 
system would have a broad spectrum of applications to urban planning. 

This study is concerned with the determination of the accuracy with 
which an urban area can be categorized, evaluation of the classifica- 
tion procedure for this purpose, and further dofinition of investiga- 
tions which should be conducted. 

The multispoctral data analysis techniques used in this study wore 
developed at the Laboratory for Applications of Remote Sensing and fol- 
lowed a standard computer training procedure utilizing manually selectod 
training fields. This preliminary study served to familiarize tho re- 
searcher with tho LARS ' s data analysis procedures and helped to formu- 
late a revised data analysis technique used in the second study, "Spec- 
trally separable training class approach to land use mapping" (Chapter 3) • 
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Expcrim ental Methodology 
Description of Test Site 

The 26.3 ha (65 acre) Meadowbrook subdivision (Figure 2 ) is located 
4 1cm. (21/2 miles) east of Lafayette, Indiana, on State Route 26 
(SE 1/4 of SW 1/4 of SE 1/4 of Sec. 23; Range Township 23N). The 
subdivision contains 68 moderately priced homos which are predominantly 
2 to 3 bedrooms, single level with attached garage. 

The subdivision is bounded on the north and west by woods and with 
cultivated fields to the south and east. Roads in the subdivision are 
asphalt with surface color from light grey to nearly black; driveways 
are generally concrete with some being asphalt or gravel. All lots havo 
good grass cover and are minimally landscaped with only a few lots having 
trees or shrub, plantings. The roofing materials are asphalt shingles, 
and their colors are predominantly light green, or blue and grey with 
several roofs being black. 

Data Collection 

The remote multispectral scanner data for this study were collected 
by the Willow Run Laboratory, University of Michigan, using a multispec- 
tral optical-mehcanical scanner mounted in a C-47 aircraft. The over- 
flight of the test site was made on May 6, 1970 at 10:30 A.M. from an 
altitude of 915 meters (3,000 feet). The data were collected in 12 
discrete channels of the spectrum (six in the visible and six in the 
reflective and emissive infrared) . A detailed description of tho scanner 
system is given by Holter ( 1967 ) . 




Figure 2. Black and white aerial photograph of the Meadowbrook 
Subdivision. 
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Data Organization and Registration 

An aIog-.Dirit.al Data c o nversion . The raw data frcra the overflight 
was in analog form and was first transformed to a digital format compa- 
tible for the LARS digital computer and multivariate pattern recognition 
programs. Each analog scan line was sampled at twice the normal digitiza- 
tion rate (440 points instead of 220) as a means of increasing the amount 
of information gained from materials on the ground. Because of a machino 
limitation bounded by the maximum rate of digitization in the analog to 
digital conversion system, doubling the number of points sampled per scan 
lino meant halving the number of channels of data from twelve to six 
(Table 1 ). Each data point, or remote sensing unit (RSU) is given a 
unique address in a two-dimensional coordinate system composod of scan 
line numbers and data points within a scan line. Each RSU is designated 
to represent a unit of area on the ground. Since these data were col- 
lected at an altitude of 915 meters (3,000 ft.), each RSU, directly 
below the aircraft, represents about 9 square meters on the ground. 

Data Registration. Because the scanner data were collected with a 

— ■ *— 1 — i ■ » — * • «»»*» — >1 ■ » ■ 1^ MMB 

three aperture scanner, one for the visible, reflective infrared, and 
emissive infrared portions of the spectrum, the respective sets of 
channels of data had to be aligned such that all data points would coin- 
cide. This process, called "Data Registration,” is accomplished by 
locating chock points on grey scale computer printouts (to bo described 
below) of the channels of data, such as road intersections or comers of 
fields. After an array of such points are established throughout the 
flightlinos, the data are aligned with those points by the computer. A 
complete description of this proceduro is given by Anuta (1970). 
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Table 1, Re-digitized roultispoctral data channels used for 
collection of data over test site on May 6, 1970. 


Channel 

Wavelength (micrometers) 

Spectral 

Description 

1 

.40 - .44 

violet 

*2 

.50 - .52 

blue-green 

3 

.55 - .58 

green 

*4 

.50 - .62 

yellow 

5 

.66 - .72 

rod 

♦6 

.80 - 1.00 

Reflective IR 


♦Channels selected by $DIVERG for the classification of the test site. 
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LARS Data Analysis 

The computer programs used at LARS to analyze multi. spectral scanner 
data were developed at LARS under a NASA grant; their functions and se- 
quential steps involved in the analysis procedure are indicated in 
Table 2 . 

Computer grey- scale printouts, which are digital, spatial displays 
of the spectral responses of the data, wore first obtained in several 
channels. From ground observation collected at tho tost site on the day 
of the flight, computer training samples were solocted from areas of 
known materials such as trees, grass, roads, driveways, and rooftops, 
located on the groy scale printouts. For this study, 14.1 percent of 
the total data points were used in computer training. The procedure 
is referred to as a "supervised" approach to computer training 
because the researcher selects only training samples of materials 
he desires the computer to recognizo for the classification. 

After training samples were selected, another program, 'COLAS, was 
used to determine spectral separability or uniformity between and among 
training classes. Often classes of materials, which wore not spectrally 
separable were merged to speed computer classification timo, and sub- 
classes of materials, which were spectrally separable, wore found that 
could greatly increase the detail of the classification. 

The next step in the analysis procedure was to obtain statistical 
information concerning the spectral responses of each class or subclass 
of materials. This was accomplished by the program, $3TA1 for which the 
output consists of class means and covariauco matrices for the responses 
in each wavelength band for each computer training class. Another program 
called $DIVERG selected those channels which ,;d tiro greatest 
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Table 2 . 


ANALYSIS 

STEPS 


1 


2 


Comnuter programs used for classification of multispec- 
tral scanner data (after West, 1972). 


PROGRAM 

NAME PROGRAM DESCRIPTION 


LARSPLAY 

$PIC Obtain grey scale printouts of each 

spectral band desired using alpha- 
numeric symbol s . Primary purpose 
to select training samples from 
known areas based on "ground obser- 
vations." 


$NSCLAS Obtain spectral groupings or clusters 

of the data based on reflectance 
values from selected spectral bands 
without supervision or prior ground 
truth designations. 


LARSYSAA 

$STAT 


$DIVERG 


^CLASSIFY 

$DISPLAY 


Statistical analysis of reflectance 
data from training fields and classes - 
histograms, spectral curves, mean 
vector, covariance and correlation 
matrices are generated. 

Selection of best subset of avail- 
able spectral bands for classifica- 
tion; obtain indication of separa- 
bility of proposed classes. 

Classification of each remote sensing 
unit in designated area. 

Display of classification by alpha- 
numeric symbols with the use of 
threshold levels. Evaluate classi- 
fication performance. 
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separation between all combinations of desired classes (Table 2 ).. 
V/acker (1971) found that tho accuracy of classification is not greatly 
increased boyond the use of four to five channels for agricultural 
materials. 

Next, tho training field statistics were used as classification 
criteria by a classifier program called ^CLASSIFY. Tho computer 
classifies each data point or RSU based on maximum likelihood criteria 
(Landgrobe, et al. , 1968), 

The final step is to display tho classified area, using tho pro- 
gram, $DISPLAY. Any alphanumeric character can be used to display a 
particular class. This program also offors an option for testing the 
accuracy of tho classification. The usor inputs tho coordinates of 
test fields for which the ground truth is known. On call, tho compu- 
ter will note tho accuracy on tho test fields and calculate the per- 
cent accuracy. 

Res ul ts _pnd Discussion 

Tho final classification of the subdivision (Figure 3) is dis- 
played in two categories: man-made, containing the classes of roof 
and streot-drivoways (Figure l '<) , and naturally occurring materials 
containing the classes troos-shrubs and grass (Figure 5)* The claos 
"s treets-driveways" in tho category of man-made materials is further 
subdivided into asphalt and concrete materials (Figure 6). 

Tabular comparisons between the computer classification rosult3 
and an estimate of the actual area covered by each category (Table 3 ) 
were obtained by overlaying a photograph of tho test area with a dot 
grid. The dots in each category were counted and dividod by tho total 





Roof Surfaces 


Grass 

Trees -Shrubs 


Streefs - Driveways 


Figure 3. 


Computer classif ication display of man-made and natural 
features of the test site. 


• 1 <i 


Roof Surfaces 


Streets - Driveways 


Figure 4. Computer* classification display of the man-made materials 
present in the test site. 





Grass 


Trees - Shrubs 


Figure 5. Computer classification display of the natural materials 
present in the test site. 


Computer classification display of the category "streets- 
driveways” subdivided into asphalt and concrete materials 
present in the test site. 
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Table 3 . 

Computer classification results compred to a dot-grid 
estimate of the different categories of material present 
in the test site (all numbers represent percent in the 
total area) . 



Man-Mado 


Material 


Computer 

Results 

Dot-grid 

Estimate 

Roof 


5.1 

8.2 

Streets- 

Driveways 


17.3 

13.3 

Subtotal 






22.4 

21.5 



Naturally Occurring 


Trees and 
Shrubs 


15.2 

20.6 

Grass 


62.1 

57.9 

Subtotal 


77.3 

7$. 5 

Grand Total 

99.7 

100.0 
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dots in the scene. The computer classification results were auto- 
matically calculated by the computer. A total of 99.7 percent of 
the data points in the test area (Table 3) are actually classified 
by the computer. The remainder of the data points were not class- 
ified because an insufficient number of computer training classes 
were present to adequately represent these data points. These points 
are considered "thresholded" points. 

Tho classification results (Table 3 ) reflect the percent of 
the total area occupied by a mapped category. The results do not 
indicate the spatial accuracy of the classification. Comparison of 
the photography (Figure 2 ) to the classification displays, shows 
that the spatial accuracy is good and that most of the misclassific- 
tion occurs near class boundaries. 

The computer classification results align themselves will with 
tho estimated actual percentages of each category as shown 
in Table 3 . The results show that the computer classification 
slightly over-estimates the percent area occupied by the category 
''man-made" (22.4 versus 21.5 percent of the area as calculated by 
the actual) and only slightly under-estimates the area occupied by 
the category "naturally occurring" material (77.3 versus 78.5 by the 
actual) . 

Within the category "man-made," Table 3 shows the percent area 
classified as roof material as 5.1- and the percent area classified 
as streets-driveways as 17.3. The estimated actual percent of those 
materials in tho total area are 8.2 for roof and 13.3 for streets and 
driveways. With closer examination of the photography (Figure 2 ), 
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concrete sidewalks and patios wore visible around many of the homes, 
but not to the extent shown by the computer classification of this 
material in Figure 6 . Much of the misclassification way' bo explained 
by shadows, and the large variation in roof size, color and surface 
orientation with respect to the scanner. 

Some difficulty was experienced in differentiating roofs from 
dark vegetation in some of the preliminary classifications, but was 
overcome by introducing a reflective infrared channel (0.8 to 1.0 
micrometers) in the classification procedure. This channel was quite 
valuable in mapping man-made materials from naturally occurring be- 
cause the responses were largely independent of color. 

Within the category "naturally occurring," the computer estimato 
of the percent of the total areas occupied by trees and shrubs and by 
grass were 15.2 and 62.1, respectively. The actual percentages of 
the total area for these two classes of material wore 20.6 and 57.9, 
respectively. Comparison of the photography (Figure 2 ) with the 
classification display (Figure 5 ) indicated that much of the under- 
estimation of trees and shrubs and over-estimation of grass may be 
duo to difficulty in discriminating between shrubs and grass materials. 

Considerable spectral variability existed within both the natural 
and man-made matorials; therefore, a large number of samples wore re- 
quired to represent single classes of materials. This is especially 
true in categories dealing with materials of small area. A possible 
explanation of the variability in these smaller areas Ls in the regis- 
tration of the data channels. By ovorlaying each individual channel 
of data (in the form of grey scale printouts of the data) on a light 
table, it is found that registration is generally within one or two 


data points betweon channels. For large area analysis, this degreo 
of registration accuracy would not generally present any problems, 
but for small areas this degreo of error can bo critical. Consider- 
ing the ground area occupied by some of the materials such as roof- 
tops, streets and driveways, it is very important that each data 
point in each data channel characterize the same exact spot on the 
ground. 

A possible improvement in classification accuracies especially 
with materials of limited ground area could be obtained if tho num- 
ber of data points for training could be increased. This could be 
accomplished by either collecting the data at a lower altitude or 
increasing tho resolution of the scanner. If the data were collected 
at a lower altitude, than tho area represented by a single data point 
could decrease. This would decrease the effective coverage of the 
scanner but would also incroaso tho resolution of the data. It 
would be most advantageous to mechanically increase tho resolution 
of tho scanner through incroased scanner technology and increase the 
effective coverage of the scanner. This would allow improved scan- 
ner resolution for more detailed studies and increase the ground 
area over which data could be collected, thus increasing tho economy 
of the scanner system. 


Concl usions 

Much morn work must bo done with multispectral analysis of urban 
land use before exact procedures and spectral bands most useful to 
this typo of analysis can bo made. However , from this limited study 
the following conclusions can bo drawn. 
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(1) Computer analysis of multi spectral data has shown to be 
accurate in characterizing and estimating a limited number of land 
use categories in a small urban scene when compared to conventional- 
ly accepted methods such as photointerpretation. 

(2) The use of the reflective infrared channels was very 
valuable in differentiating tho man-made materials from naturally 
occurring materials. 

(3) More accurate data registration would be recommended 
for this typo of analysis. 

(4) Bettor resolution of the scanner would bo dosirablo. 
Further studies using present state-of-the-arts multispectral 

scanner and computer analysis techniques should involve investigation 
into the following: 

(1) The reflective and emissive portions of the spectrum 
in an effort to further classify tan-made materials. 

(2) A more complex urban mo for the benefit of further 
testing of the multispoctral approach to urban information 
gathering. 

(3) Studies into urban-suburban relationships over time. 


CHAPTER 3 


SPECTRALLY SEPARABLE TRAINING CLASS APPROACH TO 
LAND USE MAPPING 

Introduction 

Every year land use changes occur on millions of acres in the 
United States. The greatest portion of these changes occur when agri- 
cultural and forest lands are diverted to housing, industry, highways, 
public buildings, and parks. The effects associated with these changes 
are numerous and far reaching. 

To the present, most interpretation, monitoring, and land use 
evaluation has involved the utilization of aerial photography. Because 
of its utilitarian value as a base for most resource surveys, aerial 
photography will continue to be the basic tool available to the planner, 
but as the scope of planning activities continue to expand, new means of 
data collection, handling and interpretation must become available. 
Computer analysis of multispectral scanner data is one such new technique. 

This study is concerned with determining the ability of automatic 
computer techniques in employing multispectral scanner data to obtain 
land use information. Emphasis is placed on the standardization of a 
procedure for the analysis of data with acc...ate repoatabili ty of 
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classification results. Three overflight dates of data are classified 
and evaluated in an effort to (l) determine the accuracy of the classi- 
fication results as a predictor of land use, (2) find a suitable land 
use classification scheme for use with multispectral data and (3) deter- 
mine the best date or combination of dates that will return the largest 
and most accurate land use information. 

Experimental Methodology 
Selection of Flightline for Analysis 

The data selected for this study was collected during the 1971 
Corn Blight Watch Experiment over the intensive study area located in 
western Indiana. The intensive study area was composed of thirty north- 
south oriented flightlines that extended from the northwest to the 
southwest portions of the state. Each flightline was 1.6 km. (l mile) 
wide and ranged 12.7 to 19.2 km. (8 to 12 miles) in length. 

Because multispectral scanner data were collected over each flight- 
line at two-week intervals (when weather permitted) from the middle of 
May to late September, it offered ideal coverage of a flightline for 
this study. 

The selection of a suitable flightline was made primarily on the 
basis of the diversity of land use3 px'esent. Flightline 219 located in 
Sullivan County, Indiana, was selected because it represented four 
main land use categories: agricultural, forest, urban, and water. Other 
factors considered in its selection was its relative closeness to LABS 
facilities and the researcher's familiarity with this flightline from 
previous analysis during the Corn Blight Watch Experiment (Bauer, 1972; 
Phillips, 1972; Johonnsen and Bauer, 1972). 
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Physiograrhic Sketch of Flightline 219 . Flightlinc 219 is located 
in the northeastern portion of Sullivan County (Figure 7); it is 1.6 km. 
(l mile) wide, approximately 16 kra. (10 miles) long and is oriented in 
a north-south direction. Sullivan County is located on the vest, south- 
central portion of Indiana and is hounded by Vigo County to the north, 
Greene and Clay Comities to the east , Knox County to the south and the 
Wabash River to the west. 

The topography of the north and central portions of the flightline 
is level to gently undulating with some surface drainage into small 
ponds in the northern sector. The southern portion of the flightline is 
characterized by greater relief that contains several small weat-to-east 
draining creeks. The flightline elevation ranges from 165 meters 
(550 feet) in the northern section to 138 meters ( h6o feet) in the south. 

The underlying bedrock of the area is Pennsylvanian age sandstone 
and shale of the Covemaugh series. These units are underlain by Kiss- 
issippian age sedimontaries in which extensive strip coal mine operation 
has been active to the east of the flightline (Bieberman, 19^9). 

The soils present in the flightline were formed under forest 
conditions in Illinolan and Uisconsinan age loess undorlain by Illinoian 
age glacial till. The two soil associations present in the flightline 
are the somewhat poorly drained Reesville-Iva present in the northern 
third and the southern third of the flightline and the moderately well 
drained Cincinnati -Av a- Alford present in the middle third and some of the 
southern portions of the flightline. Both soil associations have deeply 
developed profiles that are deeply leached, acid, find a lack of I 
fertility (Sullivan County Soil Survey, 1971). 



Figure 7. The location of flightline 219 in Sullivan County, Indiana 
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Climate for the area is midcontinental characterized by wide 
ranges in temperature from summer to winter. The mean daily maximum 
temperature in July is 88°}*', and the average minimum in January is 22°F. 
Average rainfall is 96.5 cm. (38 inches) annually and is well distributed 
throughout the year, but is slightly heavier in the spring and summer 
than fall and winter (Sullivan County Soil Survey, 1971). 

The soils in the flightline are used primarily for cropland, 
pasture, and woodland with cropland predominating. The poorly drained 
areas are usually in woodland- forest . The main agricultural crops are 
corn, soybeans, small grains, hay and pasture. 

Selection of Data Dates to be Analyzed 
Data from three overflight dates of flightline 219 were chosen 
from a possible eight flights made during the summer of 1971. The eight 
sets of data were collected on May 21 and at two week intervals from 
June 27 through September 28 inclusive by the University of Michigan 
aircraft. Data for analysis wore selected on the basis of data quality, 
channels of data available for analysis, altitude, field conditions, 
availability of photographic data for overflight, and time of day data 
were collected. A summary of data overflight information for the sel- 
ected dates of analysis are found in Table 4. 

All dates were first examined for uniformity of data and collection 
procedure. This included selection of those dates with the same number 
of channels of data available for analysts; those dates where the data 
were collected near the same time of day, at same altitude and same air- 
craft heading; and those dates with the same data preprocessing (i.e. , 
analog to digital processing) procedure. Those selected dates were then 
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examined fox* overall data quality and for the availability of photo- 
graphic coverage for that date. This included dates with cloud free 
data, and dates with the same flightline coverage. Finally, dates were 
selected that provided the greatest range in field conditions for the 
land use classes and were most evenly spaced over the summer data col- 
lection period. The three dates selected and meeting the above criteria 
were: May 21, July 21 and September 15 (Missions 38, 1*1 and 45, res- 
pectively) . 

Description of Flightline for the Three Data Collection Dates 

Color infrared photography (CIR) was used in Figures 8a ,b and c for 
the flightline description because it was felt that it could supply more 
information than conventional black and white or color photography . 

Also it was used because complete flightlinc coverage wa3 available for 
the dates needed, and its collection dates came closest to all the dates 
of the actual scanner data collection than any other photographic data 
available. Although the dates of the CIR photography do not coin- 
cide exactly with the actual scanner dates, the photographs provide basic 
flightline information at near flight conditions. To interpret CIR 
photography, one must remember that these materials appear near these 
colors on the CIR: vegetation , red; bare soils, green; and water, blue to 
black. These figures are summarized as follows. 

May ( Fig ure 8a ). Agricultural cropland is represented in the 
flightline for this date as bare soil (plowed) and last year's crop 
residue (sometimes weedy); it is seen on the CIR as green, l/intor 
wheat is the only agri cultural • rop that is present in the flightline 
for this date and is red on the CiR photography. Pasture and hay appear 
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Figure 8 . Color IR photographs of Flightline 219 for the three 
daba collection dates. (Scale 1:75 ,000) 

Legend 


Photograph 

Date 

Flown 

Roll 

Number 

Frame 

Number 

Mission 

Number 

a 

5/13/71 

99KF71000U 

7907 

165 

b 

7/12/71 

99HF7100U5 

8432 

175 

c 

9/21/71 

99RF710097 

103 

180 
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(a) May 


(b) July 


(c) September 



Figure 8. 





to be at the lush green stage (very bright red on the CIR photograph) and 
appear very similar in color to winter wheat. Trees and bushes are also 
red but are not completely leafed out; areas around the lake and other 
low wet areas are lined with vegetation. The urban areas where trees are 
present are most clearly visible because trees have not yet hidden many 
of the streets, sidewalks, driveways, and roofs. Highway and railroad 
right-of-ways are lined with vegetation and are clearly visible. 

July (Figure 8b ). All agricultural cropland has been planted by 
this date. Corn fields are uniformly red in the photograph and are not 
yet tasseled. Soybean fields show early and late planting dates: 
relatively young soybeans characterised by very low ground cover, and 
fields of older beans with heavy ground cover which are very difficult 
to differentiate at this stage from corn. Winter wheat fields are 
either mature (brown color on CIR) c ’harvested. Some harvested winter 
wheat fields have been apparently re . ’anted with a short cover cron and 
appear redish brown on the CIR photo. Individual roofs and other hard 
surfaced materials are hard to see in the urban areas where there is 
now heavy tree cover. 

September ( Figure 8c ). Agricultural land for this date is char- 
acterized by mostly mature crops and /ogetation. Only a few corn and 
soybean fields have been harvested. The only fields still red on the 
photograph are late planted soybeans , pasture and hay. Tree vegetation 
is still quite red especially near the ed^es of the lake. Near the 
center of the photography are several clouds and cloud shadows. (Those 
clouds and shadows appear only on the photography and not in the scanner 
data. ) 
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Collection of Ground Observations 
Black and white base line photography at a scale of 1:20,000 was 
collected along tho flightline i n April, 1971 for use in the location and 
coding of fields in the flightline. Ground conditions information was 
then collected by Agricultural Stablizations and Conservation Service 
personnel from Sullivan County A.S.C.S. office on the ground for the 
entire flightline near the time of the flight. These data included field 
information concerning location, size, identification o' cropland, 
pasturelond, forestland, illoland (land :ot presently gaged in r. - 
duction) and non- farm areas (urban, water areas, trans- -tation, ann 
extraction). This information was made available on computer printouts 
and was updated (for randomly selected corn fields) on a biweekly basis 
during the period from June 15 to September 30, 1971- 

High altitude (l8,000 meters or 60,900 feet) CIR photography 
(shown in Figures 8a, b, and c) flcwn on dates within two weeks of the 
scanner data collection was also available for general ground observation 
information. These photographs, used in conjunction with a "VARISCAH," a 
film projection device that provided the viewer with four levels of film 
magnification (3.00X, 5-00X, 12.10X, and 23.68X) , were extremely valuable 
in verification of field boundaries and ground observation information and 
in various other phases of the data analysis. 

Black and white photography collected at 5,000 feet was available 
for most scanner flights used in the analysis. This photography allowed a 
field by field cross-check of the scanner data for the presence of clouds 
and offered supplementary point data information for the ground observation 
data. It was especially useful in evaluating the different scanner data 
flights for data quality during the selection of data to be used in the 
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Selection of a Suitable Land Use Classification Scheme 

Remote multispectral scanner data relates the physical scene only 
as a function of the reflectance and emittance properties of the scene. 

The data and the LARS programs relate in no way to the spatial character- 
istics of the scene, in other words, to the size and shape of objects, 
as does photography. Therefore, a land use classification scheme could 
only be based on land uses that have spectrally unique characteristics by 
virtiie of a particular vegetation cover or material associated with a 
land use. 

The land use classification scheme was selected for use with the 
collected multispectral scanner data. The selection of the scheme was 
accomplished by incorporating sections from the classification schemes 
of the Association of American Geographers and the Canadian Land Inventory 
with additional categories not specifically outlined in either of the four 
schemes reviewed (see Chapter 1, pages 9 to 13 ). All lend use categories 
were selected on the basis of their presence within the flightline. The 
detail of the land use scheme was determined by the resolving power of the 
scanner, in other words, its ability, to "see" that category on the ground. 
As previously determined (Chapter l), the resolution of the scanner at 
1,600 meters (5,000 feet) is 5 meters by 5 meters square. A material with 
an area smaller than 25 square meters that falls within a data point can- 
not be purely characterized spectrally in a single data point. For example, 
a small road would not contribute enough spectral information to charac- 
terize a single resolution element as "road." It was found that the 
Standard Land Use Coding Manual and the Mew York State schemes were usually 
too detailed in their first and second digit levels to be compatible with 
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the collected data (see Appendix, Table B). Those categories that were 
not too detailed were also found in the two more general schemes, the 
Association of American Geographers and the Canadian Land Inventory (see 
Appendix^, Table B). 

The final land use scheme then, reflects those categories of land 
use, in theory, that can be "seen" in the data based on the resolution 
of the apparatus. Table 5 shows the selected ]and use scheme with the 
associated land use code, the ground observation information and the 
materials (components) that would comprise the land use category. The 
derived land use scheme has four first-order categories, nine second- 
order categories, six third-order categories and four fourth-order cate- 
gories. Like the four schemes reviewed, other second, third, fourth and 
fifth-order categories can be added to the scheme if needed. Point data, 
as in the New York State scheme, can also be implemented into the selected 
s cheme . 

Analysis of Data Procedure 

The previously described "training field" or "supervised" computer 
training procedure (see Chapter 2, page 19 ) has several limitations when 
more than one set of data are evaluated fo,r spectral information. 

c 

Standardizing the training procedure from one set of data to another is 
difficult because the procedure is very subjective. The computer training 
phase continues until the desired classification results have been ob- 
tained, or the "best" classification of the data has been reached. The 
training phase of the analysis is usually the more laborious; therefore, 
the time aspect is a second consideration. 


42 


• 

co 


-p 


G 


d 


C 


o 

& 


B 

O 

CO 

o 

£ 

d 

0 

G 

s 

d 

o 

o 


co 

§ 


o 

fd 

o 

§ 

w 


tz 

to 

W 

G 

o 

O 

CO 

•H 


■P 


d 

> 


P 

a) 


CO 


.o 


o 


T) 

G 


3 

o 

CO 

p 

s 

bO 

o 


H 

>d 

EH 

a; 

<5 

p 

> 

a) 

« 

rH 

w 

d 

CO 

P 

PQ 


o 

■G 



o 

d 

B 

0) 

o 

B 

(X 

d 

o 

x: 


o 


CO 


G 


O 


•H 


■P 

w 

d 

o 

o 

o 

•p 

o 

«h 


•H 

w 

CO 

CO 

a 


tu 

H 

o 

O 


d 

a 


p 

'd 

G 


3 

iH 


T) 


d 


-p 


O 


a) 


rH 


0) 

o 

co 

M 


Eh 


«* 

• 

o 

in 

M M 
co tx 


PHH 

0 

rH 

Q § 0 

■8 


E-t 

H O CO 










-p 








o 










p 



•* 





•H 










o 

#» 


rH 





P 










J? 

-p 


rH 





Pj 

V) 



e 




f. 


to 

p 


d 


"G 



P 

oS 



p 




u 


G 

o 

C 

+H 

G 

d 

G 


0 

0 



o 




o 


- o 

.C 

o 


O 

X 

O 


esc 

Jt 



Vn 




<M 


P *rt 

CO 

•rH 

fz 

•H 

•rl 

•rH 


0 

d 

cn 


•H 


w 


•H 


P -P 


-p 

U 

-p 

B 

-P 


> 


0 


C 


0 


Cl 


O d 

* 

d 

o 

d 


d 






3 




p 


<M +J 

B 

■p 


-P 

B 

43 


c 

0 




CO 

r d 



to 

•H d 

P 

d 

•H 

d 

P 

0 

c 

0 


*rl 

rH 


o 

•H 

rH 


<u 

G t? 

o 

br 

c 

W 

o 


d 

0 

Th 

to 

•H 

o 

Pc 

(A 

•H 

o 

p. 

3 d 

Vi 

d 

P 

d 


0 

r d 

u 

o 

d 

o 

•H 

>> 

0 

O 

•H 

P. 

> 

•H 

> 


> 

•H 

> 

P 

t r , 

p« 

p 

(0 


p 

*H 

CO 


■p 

s>. 

c 




G 


G 





•H 




•H 


H G 

P 

c 

H 

G 

3 

P 


Tj 

rH 

p. 

0 

CJ 

p. 

P’ 

0 

o 

p- 

-p d 


d 

-P 

d 


0 

P 

0 

rH 

o 

Jm 

0 

o 

O 

h 

0 

o 

to d 

G 

d 

to 

d 

G 

0 

d 

X 

oj 

p 

d 

P- 

p 


d 

p* 

p 

O P 

O 

P 

O 

P 

O 


> 

•r 1 


o 

PQ 

CO 

a 

o 

no 

CO 

o 

Jr; w 

ts; 

br 

K ' 

to 

^ ' 

Ui 

6 


CO 


i i 


i i i 


T 1 


d 

d 


I I I 


I I 


(C 

•§ 

p 

X 

to 


c 






g 


to 




0 




d 




to 


d 


p. 

> 



W 

•* 


>. 

rH 




C 

B 

rH 


0 



g 

c 

to to 


p 


c 



d 

3 


g 

P 

p 


S 

0 

d d 


p 

»C 

Cu 

to 


d 

■C 


d 

d 

d 

to 

r*H 


d to 

0 

d 

0 

rH 

p. 


X 

hc> 

0 

rH 


P> 

P- 


*w* 

h (0 

c 

a 

-P 


o 

X 

>. 

p 

-P 


It 

d c 

o 

0 


■p d 

•ii 

o' 

tn 


p 

o 

0 

o 

(0 


o 

-H 

p 


43 

P 

f! 


0 

0 

o 

a 

to 

to 

0 

0 

p 

P > 

CJ 

r5 

w 

t) u: 

p. 

d 

> 






> 


1 


4-3 

0 

d 

•H 


y< 

o 

u 




& 

O 

G 


X 

0 


X 'd 


O 

d 

rH 

o 




d 

H 

0 


d 

d 

0 

•H ft 

> d 

-p 

+> 

w 

(X 

(X 




W 


to? 


K 

Ph 

p. 

Cf) 

CO 


I I 


< 


PQ 


o o w 


(n O 


d 


a 

o 

•H 
■P 
CJ 

! % 
p -p 
P i rH 
p 

a o 

O *H 

£ b 

o «d 
to 

<u • 

« < 


c • 

o • 

•H • 
•P 

y p< 

3 O 
'G p 

o o 

p 

Pc 1* 

o 

Pc p 

o 

p • 
O a! 



• 






• 






• 

• 


• 



• 

• 


• 



. G • 

• 

• 

• 



O • 

• 

• 

• 



*rl • 

• 

• 

• 


c 

■P • 

• 

« 



•H 

O • 

• 


W 


a) 

G • 

• 





TJ • 



O 


bC 

O • 

d 

to 




P • 

P 

G 

P 


rH 

P- • 

3 

O 

•H 


H 


•p 

G 

4-3 

G 

d 

d S 

"J 

G 

P 

O 


wj d 

3 

•H 

O 

>r< 

w 

d ,q 

p. 

P 

O 

■P 


P 


-P G 

W 

O 

• 

O • 

• 

to O 

•H 

d 

rP 

<m d 

^1 

d CJ 

r cJ 

P 




P 


•P 


• 


O • 

• 

X 


CVJ 


0, H 

CJ 

f-1 




• 


• 




fO 


CJ 






C P 

•H 



c 

P 


c 




'd 

O 0 

" Ai 

#* 

'd 

o 

cj 

tS 

o 




d 

•H 

to P 

CO 

d 

«rl 

e 

c 


CO 

E-i 

ix 

ro 

s~s 



aJ 

p p 

+j 0 

-p 

c3 

P 


0 

p 




0 0) 

a; p. 

0 


0 

o 


eg 



't 

-p -p 

0) 

O 

•P 

-P 

0 

CJ 

P 

p 



O 0 

p •> 

P 

b 

0 

P 

0 



o 

bC > 

P to 

p CO 

O 

bu 

«M 

O 

to 


p 

,d 

0 

CO r-y 

to pi 

X 

a) 

P 

fi 

0 

C-J 

fe 


c 

w 

> £ 

1— i 

f — 1 

to 

> 

3 

CO 

> 


0 


O 

•> nj 

* a! 

•d 


CO 




p 

3 

*d 

C H 

co |* 

M M 

a 


'0 

c 

o 

o 

0 

H 

<u 

0 rH 

<p a) 

<M CD 

<L> 

0 

'0 

0 

0 

p 

Cm 

X 

<D 0 

O TJ 

O "0 

X 

<u 

P 

X 

o 

w 

0 

<m 

•H 

P .G 

O -H 

O -H 

•H 

P 

0 

«H 

p 

ts 

K 


to CO 

err w 

tc CO 

s 

i 

M 

m 


tc 

p 


1 

1 

1 

i 

1 


i 

y 


o 












CO 

















rH 

rH 







oj 

0 

to 






•H 

•H 

to 

0 





U 

P 

cd 

1* 





<U 

0 


1 





■P 

P 

to 

<H 



P 


cd 

0 


o 



e 


6 

e 

* 

1 



0 




to 

p 



e 


a; 

0 

£ 

.0 



p 


a 

o 

0 

bO 



0 


cd 

0 

P 

•r~i 



0 

'd 

<M 

«M 

,G 

P 



P 

d 

Jh 

P 

0 

to 


CO 

P CO 

cj 

d 

0 




E 

C 

■ — 1 

to 

CO 

#» 

0 0 

CO to 

0 

0 O 

. 



to 

I* > 

<U '0 

4) 

P O 

to 

*d 

'C 

<D 

PI -C 

>- c 

P 

0 bC 

Jh 

Jh 

P 

a> 

bn tc 

0 O 

P 

0 0 

cj 

S 

0 

Jh 

•rH »H 

t-0 Pu 

CO 

5= P5 

* “ 

w 

« 

e- 

« w 


§ 

u 

W M 
CO 

m a 

r A 5T 


^ O CO| 



• 

• 

• 

• 


• 

• 

• 

• 

• 




• 


• 

• 



• 

p 


• 

• 

• 

• 



• 

c 

• 

• 

• 



• 

0 

• 

• 

• 



0 

E 

• 

• 

• 



M P 

P 

• 

• 

• 



0 0 

0 

• 

• 




PI P 

0 

• 

• 



0 

0 

P 




G • 

o 

TJ ^ 

P 


0 


O • 

P 

0 


"d 

p 


•H 

0 

0 ,C 

0 

d 

o 

s 

p b; 

o 

W 

P 

cj 

o 

as 

0 C 

CO 

>0 0 

CO 

rH 


rO 

p "H 

0 

0 P 

0 

rd 

c 


P P 

PS 

T0 S-l 

> 

to 

0 


o o 


C 



p 


Q P 

p 

O • 

m 

cj 

C P 


M O 

0 

D-. rH 

CVJ 

Si’! 

0 P) 

CO 

C S' 

p 




PI 


0 

0 

• 


• 

P * 

• 

P • 

r» 

<t 


« 

0) < 


Eh < 

• 




a 


• 

H 




M 


> 

M 




M 


M 


i 



M . 




Railroad N -Railroad right-of-way -Mixed short and tall 

green vegetation 


Anticipating these handicaps of the "supervised" computer training 
procedure, a more objective method of training was selected. The concept 
of determining spectrally separable classes within a set of data was 
decided as the objective approach to be taken in evaluating the three 
sets of data. Ideally, if the maximum number of spectrally separable 
classes for a set of data could be determined, then the maximum amount of 
information could be obtained from a set of multispectral data. The LAR3 
program $IiSCLAS (W acker, 1969 ) was selected as the most suitable program 
to analyze a set of data for spectrally separable classes. 

The analysis of data procedure (Table 6) was identical to that 
described in the preliminary study (see Chapter 2, page 19) with the 
exception of the substitution of a "non-supervised" classifier program , 
$N3CLAS, for the "supervised" selection of training fields. The $N5CLAS 
program essentially "provides the user with the capability of classifying 
a limited number of data points on a nonsupervised basis" (Wacker and 
Landgrebe , 1971). It is nonsupervised in that the researcher is not in- 
volved in the process of selection of computer training classes. The 
user must identify the classes after clustering is completed. 

To cluster (or classify) a set of data points, the user must first 
select the desired number of classes to be found in the data. The points 
from the specified area of data are read from Aircraft Data Storage Tapes 
into the computer for processing. After the desired number of classes 
are found in the data a map displaying the clustered area can be dis- 
played on the line printer. Tables containing the means and variances of 
each class as well as the pairwise separability values between all class 
pairs are listed on the printer output. 
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Step 

1 

2 

3 

4 

5 

6 

7 


Table 6 . Analysis of data procedure, 
(with program options used) 

Program or 
Procedure 


Onsclas 


Regroup $NS 
classes 

vSTAT 
$DIVERG 
$CLASS 
$ DISPLAY 


Assignment of 
classified classes 
to land use categories 


Options 


MAX. CLASS = 35 

CONVERGENCE = 100$ 
MIN. POINTS = 2 

SEPARABILITY * 1.0 


BEST 4 CHANNELS 
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$NSCLAS Procedure 

All twelve data channels were used in the $NSCLAS an.alysis in order 
to maximize the use of all the available information collected. 3y 
knowing the program requirements for computer core (Wacker, 1969), using 
all twelve channels of data, it was calculated that about 14,200 data 
points could be sampled from the flightline. Assuming that the flightline 
is 180 columns wide 1 and 1000 lines in length, the sampling rate is equal 
to 2.5 percent. 

The program options used in the data analysis were selected for 
ease of data handling and conservation of computer time. Two of the three 
options used were "MAX. CLASS" equal to 35 classes and "C01TV" (conver- 
gence) equal to 100 percent. "MAX. CLASS" simply directs the program to 
find 35 spectral classes (or clusters) within the data. "C0NV" is an 
option for the percent of data points that must meet program requirements 
to be used in a class assignment within the program. It was felt that 
the extra information that may be gained by using the maximum number of 
classes available in the program, which is 4o classes, with the same con- 
vergence or less (100 percent or less) would not Justify the increased 
computer time that would be required. 

Because of a limitation in the LARS program OdIVERG it was required 
to reduce the number of classes to a maximum of 18 classes or less for 
classification. With the use of the class separability information pro- 
duced by the $i'ISCLAS program a procedure was established to reduce the 


1 The width of the flightline was reduced by 20 columns on each side 
because of errors that may be encountered duo to scanner "look angle" 
effect (Tanguay, l*69c). 
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initial 35 classes to 18 or less. Because $IiSCLAS has the tendency to 
"find" the desired number of data clusters (classes) in the data which 
are most easily separated, it was deemed desirable to allow the program 
to first find the maximum number of classes in the data then reduce these 
classes to 18 or less manually by using a fixed separability. 

Procedure for Recombining Clustered Classes 

Using the ordered list of class separability information from 
5NSCLAS, two classes were combined if their separability value was less 
than 1.0. The separability value of 1.0 was used because it related a 
ratio of unity for the distance between a pair of cluster (class) centers 
and the sum of their cluster spread or variance. In theory, when the 
ratio of the distance between the class centers, D (Figure 9), to the 
sum of the class variances, r + r' , is equal to 1.0, the outer limits of 
the two cluster classes are touching. For separability values between a 
pair of classes of less than 1.0 the classes were termed spectrally 
inseparable and grouped; for those pairs of classes with separability 
values greater than 1.0 the classes were termed spectrally separable and 
not grouped. 

The regrouping procedure for the $NSCLAS derived classes started 
with the pair of classes with the lowest separability valtie from the 
ordered separability list and stopped with that pair of classes with a 
separability value of 1.0. A pair of classes were not grouped together 
until all member (linked) classes were linked to each other. Figure 10 
shows a simplified regrouping procedure for only six classes reduced to 
three. Classes A, B, C and E, F were grouped because each member class 
was linked with each other. Group E, F was not grouped with group A, B, C 



sponse in Channel Y 





4 9 



Figure 10. 


R e _,T r0 uping procedure for reducing six classes to three, 
(numbers represent separability values between classes) 
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because not all member classes were linked; for the same reason class D 
cannot be grouped with A, B, C. 

This procedure was applied to each of the three dates of data. 

The 35 $NSCLAS generated classes for the May, July and September data 
dates were reduced by the above outlined procedure to 12, 16 , and 12 
classes, rsspectively. Th8 mere fact that the May and September data 
reduced to a fewer number of classes than did the May indicated that 
the May and .September data contained fewer spectrally separable classes 
than did the July data. This was expected because the May data were 
characterized by mostly bare soil, some green vegetation and water, 
and the September by mostly mature or harvested crops and water; how- 
ever, the July data reflected many complicated green vegetation 
differences, bare soils, and water. 

Tlie third option used, MIN. POINTS, equal to 2, provided a means 
by which the coordinates (line and column number) of two or more data 
points that belonged to the same yi'TSCLAS class were saved on computer 
punch cards for input into the next computer program in the analysis , 
$STAT. The MIN. POINTS option of 2 or more data points was used for 
selection of every data point because of the high variances usually 
associated with classes composed of single data points. Errors in data 
overlay (registration) and lack of nurity of single data points explained 
most of the high variances. 
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Results and Discussion 

Interpretation of classification results is accomplished by 
randomly selecting representative test areas of different land uses from 
the multispectraJL data obtained for each flightline date. From these 
test areas the total number of land use sample points in each different 
spectrally separable computer class were tabulated. This tabulation for 
each flightline date yielded the percent correct recognition of a land 
use category by different spectrally separable computer classes (Appendix, 
Tables Bl, B2 , and B3). 

The assignment of a land use category to a spectrally separable 
computer class is made by selecting the computer class or classes with 
the highest percentage of samples identified with that land use category. 
In some cases , a land use category was not associated with any specific 
computer class and therefore, could not be differentiated for that flight 
date. Conversely, several spectrally separable computer classes often 
mapped the same land use category. In this instance the computer classes 
were grouped and their results summarized for the specific land use cate- 
gory. A summary of the grouped spectrally separable classes for the 
assigned land use categroies for their flight dates is found in Table 7. 

Only one land use category was assigned to a spectrally separable 
computer class because only one alphanumeric display symbol can be used 
to display each computer class. This symbol was assigned to each computer 
class to graphically represent its distribution on a computer printout 
display for each classified flightline date (Figures 11a, b, and c). For 
grouped computer classes (i.e., computer classes that are assigned to the 
same land use category) the same symbol is assigned to each computer class 
within the group. 
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Table 7. Summary of land use category assignments to $HSCLAS derived 
spectrally separable computer classes. 

Land 

Use 

Land Use Code Spectral classes 1 


May .Tuly Sept . 


I. Resource Production 

A. Agricultural 

1. crop production 

a. row crop A 

b. small grain . . . . B 

2. forage production 


a. hay C 

b. pasture D 

B. Forest 

1. continuous ...... E 

2. discontinuous . . . . F 

C. Extraction G 

II. Water Resource 

A. Ponded and Lake 

1. fresh water H 

2. waste treatment ... I 

B. Marshland J 


2,3,4 

6,11 

3,7 

7 

2, 3 

1,4 

9 

7 

5,6 

10 

12 

— 

11 

9,13 

10,11 

— 

l4 

— 

— 

4, T* 

T 


8 

16 

— 

12 

15 

12 

— 

8,10 

8, : 


III. Urban 

A. Urban Core K 

B. Suburban L 


5,6, T 1 2 


1 Numbers represent the unique spectrally separable classes for each 
flight date. 

* The letter "T" represents data points that were "thresholded" 

(page 25). 


Figure 11. Computer land use classification displays and actual land use categories present 
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By computing the percentage of each laud use category that is 
associated with each computer class or grouped computer classes , an in- 
dication of the recognition of each land use class is obtained (Tables 
8a, 8b, and 8c). Land use categories that are not recognized by any 
computer class appear as columns of dashes. Ideally a well-mapped land 
use category is characterized by none to very few low percentages of 
"other" land use categories appearing in the same column. An example of 
a well-mapped land use category is land use "I" (waste treatment) in the 
May and July data flights (Tables 8a, and 8b, respectively). Although 
the accuracy of land use category "I" for the July flight is much better 
than that of the May (91% vs. 59/0, both dates show that no other land 
use categories are confused with the recognition of land use "I". 

The accuracy of the spectrally separable computer classes in re- 
cognizing a land use category is calculated by totaling the percentages 
of each assigned land use category in each computer class for each flight 
date. This information plus the average recognition accuracy by land use 
category and by flight date is summarized in Table 9; a summary of the 
"best" computer recognition results by land use category considering all 
three flight dates is shown in Table 10. 

Classification Results by Land Use Category 

Agriculture: crop production: "row crop" (Code A ). The crops corn 
and soybeans represent the land use category "row crop." It is mapped 
best in the May flight (Figure 11a) at an accuracy of 99 percent. This 
accuracy can be attributed to the fact that all the land representing 
this category is predominantly bare soil for this date. Kristof and 


able 8. Percent recognition of land use categories with multispectral techniques. 
< a ) May (b) July 





sofuoSoquo osfl puuq xunqov 



ee Table 7 for definition of land use categories. 
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Table 9. Summary of 
results. 

percent recognition 

Of i 

computer classification 

Land Use 
Category 

Land Use 
Code 

May 

July 

Sept. 

Row crop 

A 

99 

25 

49 

Small grain 

B 

23 

6o 

62 

Hay- 

C 

38 

99 

73 

Pasture 

D 

39 

11 

— 

Forest , continuous 

E 

95 

56 

94 

Forest, discontinuous 

F 

— 

90 

— 

Extraction 

G 

— 

77 

94 

Fresh water 

H 

4o 

2k 

— 

Waste water 

I 

95 

91 

100 

Marshland 

J 

— 

6o 

75 

Urban core 

K 

68 

ko 

25 

Suburban 

L 

— 

9 

— 

Average for Mapped Categories 

62 % 

5k% 

11% 

Overall Average 


hl% 

5k% 

kl% - 
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Table 10. Best overall computer classification results using all dates. 


Land Use 
Category 

Land Use 
Code 

Date 

Percent 

Accuracy 

Row crop 

A 

May 

99% 

Small grain 

B 

Sept. 

62% 

Hay 

C 

July 

99% 

Pasture 

D 

May 

39% 

Forest , continuous 

E 

May 

99% 

Forest, discontinuous 

F 

July 

90% 

Extraction 

G ' 

Sept . 

9b% 

Fresh water 

H 

May 

ho% 

Waste water 

I 

Sept . 

100% 

Marshland 

J 

Sept. 

79% 

Urban core 

K 

May 

68% 

Suburban 

L 

July 

9% 


Overall Average 


12 % 
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Zachary (1971) have shown that computer analysis techniques using multi- 
spectral data have mapped bare soil from other surface covers. How crops 
are the only land use which are found as bare soil at this time of year 
in Indiana. 

For the other two dates row crops were not classified very well. 

The prime reason appears to be that this land use category is not char- 
acterized by crop types which are clearly distinct from other categories 
(Tables 8b and 8c). The literature reported that remote sensing 
techniques can distinguish well between these crop types (Chapter 1, 
page 5) . These results were arrived at through thorough computer 
training on specific crops. The techniques utilized in this study 
did not allow such precise computer training. As this procedure is 
refined, higher recognition accuracies between crop types would be 
expected. 

A griculture: crop production: "small grain" (Code B ). This land 
use category is represented by ten fields of the cover type, winter wheat. 
The best results (60 and 62 percent for July and September, respectively) 
are obtained when the crop is harvested and the fields are in stubble. 

The July photography (Figure 8b) shows indications that some of the winter 
wheat fields have been replanted with a cover crop which is still green 
in September (Figure 8c). Apparently the mixture of stubble and light 
cover croj) for the July data and the abundance of a green cover crop in 
the September data improve the mappable characteristics of this class for 
these dates. This information helps to explain that wheat is confused 
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with one other category for each date: "hay" in July and "row crop" in 
September (Tables 8b and 8c, 32 and 17 percent, respectively). 

The poorest identification accuracy for this class (23 percent) 
was obtained from the Kay data. The greatest confusion associated 
with the class was "row crop" (26 percent). This confusion is probably 
explained by several areas of low ground cover in many of the wheat fields 
for this date, and these areas are mapped "bare soil" as is the category 
"row crop" for this date. These is also some confusion with other cate- 
gories which are associated with other lush, green cover materials, in 
particular, the categories "hay" and "pasture." Another possible ex- 
planation to this confusion is that the ten small wheat fields did not 
contribute enough training samples to adequately train the computer. For 
most of the other categories much larger fields were available for train- 
ing. 

Agriculture: forage production: "hay" (Code C ) . Land use category 
"hay" is mapped best (99 percent. Table Ob) in July and next best (73 per- 
cent, Table 8c) in September. The July classification is mapped well and 
not confused with any other land use category, but the September data are 
confused (26 percent) with the category "marshland. " In reviewing the 
photography (Figure 8c) much of the short, green marshland vegetation 
appears to resemble many hay fields. The May data yields the poorest 
accuracy (38 percent) because of much confusion (53 percent) with the 
dense forest class: "forest, continuous." Two possible explanations for 
this confusion are that both categories are quite green at this date and 
the possibility that the tree cover is not dense enough such that the 


grass around the trees is seen. 
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Agriculture: forage production: "-pasture" (Code D ) . Land use cate- 
gory "pasture" is not mapped very veil for any date. It is mapped only 
for the two dates May -and July with accuracies of 39 and 11 percent, 
respectively. In the May data "pasture" is confused with the classes 
"hay" and "forest, continuous," 30 and 13 percent, respectively. The 
July data show that the same confusion exists hut to a much greater 
degree; the category "pasture" is related to eight other categories 
(Table 8b). 

Although this land use category is not assigned to any spectral 
classes in the September data, it is still possible to observe the other 
categories in which it is misclassified. The largest amount of this mis- 
clas si fi cation occurs in the categories "hay" and "row crop." From the 
photography for this date (Figure. 8c) it is seen that ground conditions 
for pasture are very similar to some fields of corn and soybeans (class 
"row crop") and hay. 

Forest: "continuous" and "discontinuous" (Codes S and F ) . The land 
use category "forest, continuous," which is composed of heavily forested 
land, is mapped best in May and September at 95 and 9*t percent accuracy, 
respectively (Tables 8b and 8c). The classification on both dates is not 
confused greatly with any other categories. The July data map this cate- 
gory with considerably less accuracy (5 6 percent) because 4l percent of 
this category is mapped in the less densely forested class, "forest, 
discontinuous . " 

The category "forest, discontinuous" is mapped in July with an 
accuracy of 90 percent with only 10 percent assigned to the category 
"forest, continuous." For the dates May and September (Tables 8a and 8c) 
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when this class is not mapped, the majority of its points are mapped into 
the category "forest, continuous." 

"Extraction" (Code G) . The category "extraction" is composed of two 
small areas, totaling about 0.3 km. 2 , located in the center portion of the 
flightline. The class is characterized by spoil, swampwater areas, and 
low shrub and forested areas. It is mapped best in September at 94 per- 
cent accuracy. Little confusion from other categories is associated with 
this category for this date; it is displayed as "blanks" in the computer 
printout (Figure 11c). A small percentage of other materials are also 
displayed as blanks on the computer display. These data points represented 
materials from several categories such as "urban core" and "extraction" 
that were not recognized by the computer during classification. These 
data points are referred to as "thresholded" points and are represented by 
the letter "T" in Table J. 

The July data show the second best results of 77 percent for this 
category, and errors showed confusion with the categories "fresh water," 
"small grain" and "urban core" (Table 8b). Considering the overlapping of 
materials included in these categories, it is reasonable to expect con- 
fusion from these categories. 

For the May data this land use category is not mapped well enough to 
be assigned to a computer class. It is confused with several categories 
such as water, bare soils, green vegetation and hard surfaced materials. 
Spring rains have raised the water level in the areas of "extraction," and 
and computer has mapped them accurately as "water" (Figure 11a). 
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Water Resource: ponded and lake: "fresh water" and "waste treatment 
( Codes H and I ) . The category "fresh water" is consistently shown to be 
poorly mapped or not mapped at all, whereas the class "waste treatment" 
is consistently mapped well. Studying Tables 8a, 8b, and 8c, it is seen 
that there is a close relationship between these two categories. For 
each of the three dates where the major percentage of "waste treatment" 
is mapped there is an associated high percentage of "fresh water" also 
mapped: May, 100 and 99 percent; July, 98 and 71 percent, and September, 
100 and 99 percent. From these results it is suggested that the two 
categories should be grouped together for these sets of data. In the 
selected computer classification displays (Figures 11a, b and c), these 
categories are grouped. 

Water Resource: "marshland" (Code J ). Marshland is characterized 
by short and tall grasses, some trees and some water covered areas. This 
category is only mapped in two dates (July and September) with 60 and 75 
percent accuracies, respectively. The July data are most associated with 
the category "extraction" (27 percent) and the September are most asso- 
ciated with the category"forest , continuous" (18 percent) and to a lesser 
extent with the category "pasture" (7 percent). One would expect this 
confusion due to the overlapping of materials with other categories simi- 
lar to the land use category "extraction." 

Urban: "urban core" (Code K ). "Urban core" is characterized by man- 
made hard surface materials and an almost complete lack of green, natural 
materials. This class is mapped for all three dates. This category for 
the May data, mapped with an accuracy of 68 percent, is confused most with 


"row crop" and "small grain” 


those categories identified largely because 
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of the reflectance of bare soil. For July this category is mapped with 
an accuracy of. 40 percent and is confused with two categories that have 
some areas of bare soil or highly reflective surfaces such as "extraction" 
and "small grain." "Urban core" for September is mapped with only 25 per- 
cent accuracy. Most of the confusion (49 percent) comes from the cate- 
gory "small grain" with "row crop," "hay" and "extraction" categories 
accounting for the remainder 

Urban: "non core" (Code L ). The land use category (urban) "non 
core" is mapped only in the July data and there with a poor accuracy of 
9 percent. For this date this category is confused with nine other 
classes that relate to every land use category except "fresh water" and 
"waste water." Similar confusion is also found in the other two dates. 

The reason for this is that an urban scene is a very complex area to char- 
acterize; its range in materials may include all types of natural 
and man-made materials. This explains why this category is associated 
with so many different categories. 

Transportation: "motoring" and "railroad" (Codes M and N ) . The 
table of results (Table 9) show that no test results are displayed for 
these two categories. The physical area of materials representing these 
categories is too small to be purely characterized in either computer 
training or test samples. The actual area covered by materials associated 
with t^ese categories is one reason why these categories could not be 
mapped; the other is a problem in the registration of data channels. The 
overlay of data points in each data channel (data registration) is usually 
within ± 2 data points. In other words, the location of a data point in 
a particular channel of data would be aligned within ± 2 data points in 
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another channel of data. Because of the small sized areas of these cate- 
gories and the errors associated with data registration, it is very dif- 
ficult to obtain pure samples for computer training and testing. As is 
suggested in the preliminary study conclusions (Chapter 2, page 27 ) » 
collection of data from a lower altitude or an improvement in scanner 
resolution would have to come about before good results could be obtained 
for these land use categories. 

Results of Best Classification for Combined Dates 
Using the best classification results for each land use category, 
an accuracy of 72 percent for all land use categories was obtained 
(Table 10). No combining or grouping of categories were used to arrive 
at these results. Had computer classes been grouped to map combined 
land use categories, the accuracy would have been higher. 

An alternative to using only the best single classification result 
for each land use category would be to use a temporal data overlay 
procedure. This procedure would produce data such that each data point 
is overlayed by the corresponding data point from each successive date 
of collected data. The concept behind such a procedure is that the 
information contributed by each successive data. point over time is 
additive. Analysis of the "layered" data could then possibly yield 
better results than using any single set of data. 

Evaluation of the Selected Classification Scheme 
A minimum level of accuracy in the interpretation of the data of 
about 85 or 90 percent or better is suggested by Anderson (1971b). This 
was suggested because results such as these would be nearly comparable 
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with the level of accuracy attained by the Bureau of the Census of 
Agriculture. Although, for any date analyzed, the overall accuracy did 
not reach 85 or 90 percent; the 5U percent accuracy for July (Table 9) 
is impressive considering the fact that these results were arrived at by 
automatic methods. 

It is felt that the lack of spectrally separable computer classes 
contributed in two ways to the loss in accuracy of the land use 
classification scheme. First, it left four land use categories 
unassigned for both May and September data. The classification 
accuracy of the mapped categories for these dates (Table 9)> are 62 and 
71 percent, respectively, but when the overall accuracies are 
calculated for these dates, the accuracies dropped to 4l and hj percent, 
respectively (Table 9)» Had more spectrally separable classes been 
found in the data for these dates , the four unmapped land use categories 
for each date probably would have been mapped. Secondly, it crowded 
some land use categories together and caused confusion and decreased 
classification accuracies. Increasing the number of computer classes 
would essentially increase the possibility of land use categories 
falling into single computer classes. 

Evaluation of the Spectrally Separable Training Class Approach 
To a Land Use Classification System 

The overall design of this land use classification approach is 
centered around the concept that spectrally separable computer classes 
are determined for a set of data. These classes are then interpreted 
and used to identify and map those land use categories which they 
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represent in the multispectral data. Two problems seemed to prevent this 
system from operating more effectively. The first problem was the 
physical overlapping of componets within land use categories. The 
second, whichmay be in part a function of the first problem, was the 
apparent inability for $N£CJLAS program to identify spectrally separable 
computer classes that could uniquely characterize single land use 
categories. 

Physical overlapping of land use categories seemed to be from the 
fact that components within some land use categories did overlap 
within a land use scene. For example, the categories "forest land" 
(trees) and "pasture land" (grass) are both found within the category 
"urban". Little can be done to correct misclassification within these 
categories. One possible solution is to revise the classification scheme 
to eliminate the category overlapping, but this would tend to decrease 
the detail of the present classification scheme. Another possible 
improvement in differentiating between categories would be to introduce 
spatial data with the spectral data. The techniques that are presently 
used for this type of analysis are relatively new and for the most part, 
are still in developmental stages. 

The second problem involved land use categories such as "row crop" 
and "small grain" or "forest" and "hay", that (during the growing 
season) should be spectrally separable from each other. Ideally, each 
spectrally separable computer class ($NSCLAS selected classes) should 
characterize only one land use category; in very few cases is it found 
that a land use category is identified "purely." In other words, no 
other category should be found in a column or row for a particular land 
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use category (Tables 8a, 8b, and 8c). Each of these three tables should 
have a diagonal of highly correct identification percentages for each 
"identified” land use category; it would be said then, that each land 
use category was mapped "purely." In the July and September data 
(Tables 8b and 8c), this is not the case. The overlapping or 
confusion between these categories seems to relate to the problem of 
the program $NSCLAS not being able to determine spectrally separable 
computer classes that adequately characterize individual land use 
categories in the data. A closer look at the July data provides an 
insight into some possible improvements in the procedure for analyzing 
the data. 

The July data indicate an advantage for determining a larger 
number of spectrally separable classes in the data. In July (Table 7) 
all land use categories are assigned to at least one or more spectrally 
separable computer class. For each of the other dates. May and 
September (Table 7) four land use categories are unassigned. The increase 
of four computer classes in the July data does indicate an increase in 
the amount of spectrally separable information for that date. The 
question is raised, however, as to whether the increase in the number 
of mapable land use categories is attributable to the greater number of 
computer classes or to the increase in spectral information or both. To 
answer this question with the data collected from this study would be 
difficult, but one can speculate as to the effect of increasing the 
number of computer classes. Increasing the number of computer classes 
would spread the related land use categories among more computer classes. 
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Another possible means for increasing the separability between 
land use categories (and thereby increasing classification accuracy) 
would be to increase the sampling rate of the flightline. The sampling 
rate used for this study (2.5 percent of the data present in the 
flightline) was determined by $USCLAS program limitations, (page 46). It 
is quite possible that this rate did not allow adequate sampling of 
small and heterogenous land use areas. This, in fact, partly explains 
the reason why the land use categories within "Transportation" could not 
be mapped (pages 64 and 65 ). 

In summary, effects of some of the problems associated with the 
spectrally separable computer class approach to land use mapping can be 
possibly minimized in the following ways: 

1. increase the number of spectrally separable classes from 
$USCLAS. 

2. increase the sampling rate of the flightline. 

3. revise the classification scheme. 

4. include spatial data with the spectral data in land use 
analysis procedure. 

A Proposed Land Use Classification System 
As already mentioned, the results seem to indicate confusion in the 
classification system from two areas: one, the program $NSCLAS and, two, 
the classification scheme itself. Two changes have already been proposed: 
to increase the number of spectrally separable classes from $IJSCLAS and 
increasing the sampling rate of the flightline data. A second change 
will involve the rearrangeing of the analysis procedure and of the 
classification scheme. The change in the analysis procedure is discussed 


first. 
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The confusion between categories within the selected classification 
scheme has been shown to be caused by the overlapping of materials or 
even complete land use categories. It is suggested that the step in the 
analysis procedure where land use categories are assigned to spectrally 
separable computer classes (Table 6) be altered to assign ’’basic land 
characteristic materials" to the spectrally separable computer classes. 
The idea behind this change is that one or more "basic materials," such 
as grass , trees, bare soil, water, hard surface materials, and etc., are 
present in varying percentages in the make-up of any land .use scene. 

The object is to spectrally characterize these classes of materials by 
assigning them to the $NGCLAS determined classes. Once this step is 
accomplished a new classification scheme must be derived. 

The results seem to point toward a "ratio" approach for automatic 
land use classification. Many of the land use categories are not "purely 
mapped by their assigned spectrally separable computer class , but are 
often mapped by percentages of many land use categories. 

By assuming that each land use category can be characterized by a 
particular ratio of "basic material" classes (grass, trees, bare soil, 
etc.), a program could be written to classify a set of data, grid cell 
by grid cell, 1 into appropriate land use categories of similar ratios of 
"basic materials." To classify, the program would tabulate the percen- 
and ratios of the "basic materials" present in each grid cell and match 
that cell with the land use category that would have the closest ratio 

1 A grid cell would simply be a single unit within a set of gridded data 
that would be a given number of data columns wide and a number lines 
long. The cell size could vary as the detail of the classification 
varied. 
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of "basic material." For instance, the land use category "urban" may 
be characterized by ranges in ratios of materials of 1+0 to 90 percent 
"hard surface" materials, 0 to 30 percent "trees" and 0 to 30 percent 
"grass." Should the calculated ratios of these materials for a parti- 
cular grid cell fall within the designated range of the above given land 
use category "urban," then that grid cell would be classified in the 
land use category, "urban." The designated ranges of materials for a 
given land use category would vary with the time of year. For example, 
land use category "agriculture" could vary from ratios of "bare soil" 
"grass" and "trees" in the spring to no "bare soil" materials and ratios 
of specific "crops," "grass" and "trees" for the summer. 

The flexibility of such a system as this is a unique attribute. 
Although the analysis procedure of the proposed system is patterned after 
a land use classification study, it is applicable to a wide selection of 
users. The interpretation of the $NSCLAS spectrally separable classes by 
the user (Stage 3, Table ll) for the basic information in which he is 
interested is what determines the use of the data. 

This system would allow preanalyzed data (Stage 1, Table 11 ) to be 
available and usable to a large cross-section of users. The user would 
complete stages 3 and 5 (Table 11 ) for his own specific use. This system 

could save valuable computer processing and analysis time presently needed 

> 

for multispectral data analysis. 
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Table 11. Analysis steps in the proposed new land classification system. 
STAGE DESCRIPTION 

1. Collection of Multispectral Scanner Data 

2. oNSCLAS Processing - to determine spectrally separable 

classes within data. 

3. $NSCLAS Interpretation - assignment of computer classes 

to represent "basic" scene materials. 

4. Classification of Data - based on spectrally separable 

classes as computer training fields. 


Specific User Program - (Land Use Analysis Program) it 
will ratio the identified "basic" scene materials 
in a data cell and assign it to the closest land 
use category of the same ratio. 


5 . 
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Conclusions 

Although much work is yet to be completed before a completely auto- 
mated land use classification system can hope to be operational, the 
performance of this first attempted land use classification system uti- 
lizing only multispectral scanner data to characterize land use categories 
has produced these conclusions: 

(1) The results indicated some limited ability for automatic 
computer analysis of remotely sensed scanner data to characterize 
broad land uses within an agricultural scene. 

(2) The results indicated a need for a change in the data 
analysis procedure and land use classification scheme. 

Further study is suggested in the following areas of multispectral 
scanner data analysis: 

(1) Proceed with the development of a "ratio technique" for 
land use classification. 

(2) Further the development of the "unsupervised, spectrally 
separable class approach" to computer training utilizing the program 
$NSCLAS , 
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APPENDIX A 


Table A . A COMPARISON OF SELECTED LAND USE CLASSIFICATION SCHEMES 

(after -the Conference on Land Use Information and Classifi- 
cation, 1971). 


Scheme from Study by Association 
of American Geographers 

I. Resource Production and Extrac- 
tion 

A. Agricultural 

1. Crop production (cropland) 

2. Fruit (orchards, groves, 

and vineyards) 

B. Grazing 

1. Rangeland grazing (range- 
land) 

2. Livestock pasturing 
(pasture) 

C. Forestry 

1. Commercial 

2. Non- commercial 

D. Mining 

E. Quarrying 

II. Transportation, Communication, 
and Utilities 

A. Transportation 

1. Motoring (highways, parking, 
terminals, etc.) 

2. Railroading (rights-of-way, 
yards, terminals, etc.) 

3. Flying (airports) 

4. Shipping (inland waterways 
and marine docks and related 
facilities) 

B. Communications 
1.. Telephone lines and facilities 

2. Telegraph lines and facilities 

3. Radio stations and facilities 

4. Television stations and 
facilities 


Canada Land Inventory 
I. Urban 

1. Built-up area 

2. Mines, quarries, sand and 
gravel pits 

3. Outdoor recreation 

II. Agricultural Lands 

1. Horticultural, poultry, and 
fur operations 

2. Orchards and vineyards 

3. Cropland 

4. Improved pasture and forage 
crops 

5. Rough grazing and rangeland 

a. Areas of natural grass- 
lands, sedges, herba- 
ceous plants, and aban- 
doned farmland whether 
used for grazing or not 

b. Woodland grazing 

III. Woodland 

1. Productive woodland 

2, Non-productive woodland 

IV. Wetland (swamp, marsh or bog) 

V. Unproductive Land (land which 
does not, and will not, sup- 
port vegetation) 

1. Sand 

2. Rock and other unvegetated 
surfaces 


VI. Water 


Table A (cont.) 


C. Utilities 

1. Electric- 

2 . Gas 

3. Water (including irrigation) 

4. Sewage disposal 

5. Solid waste disposal 

III. Urban Activities 

A. Urbanized livelihood areas 
(urbanized 'areas defined by the 

Bureau of the Census) 

1. Industrial 

2. Commercial 

3. Services 

4. Residential 

5. Recreational 

B. Other urban livelihood (places 

of more than 2,500 population but 
not including urbanized areas 

1. Industrial 

2. Commercial 

3. Services 

4. Residential 

5. Recreational 

IV. Towns and Other Built-Up Livelihood 
Areas 

A. Industrial 

B. Commercial 

C. Services 

D. Residential 

E. Recreational 

V. Recreational Activities (other than 
those in urban areas and towns) 

A. Mountain oriented 

B. Water oriented 

C. Desert oriented 


D. Forest oriented 


Table A (cont.) 


E. Other (including combinations 
of above) 

VI. Low-Activity Areas 

A. Marshland oriented 

B. Tundra oriented 


C. Barren land oriented (including 
lava flows, dunes, salt flats, 
mountain peaks above timber 
line, etc.) 

VII. Water-Using Activities 

A. Lakes 

B. Reservoirs 


C. Streams 
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Table A (cont.) 


Standard Land Use Coding Manual 
(One- and Two-Digit Levels Only) 

1. Residential 

11. Household units 

12. Group quarters 
Residential hotels 
Mobile home parks or courts 
Transient .lodgings 
Other residential, 

NEC- 


13. 

14. 

15. 
19. 


2. Manufacturing 


21 


Food and kindred products - 
manufacturing 

22. Textile mill products - 
manufacturing 

23. Apparel and other finished 
products made from fabrics, 
leather, and similar materials - 
manufacturing 

24. Lumber and wood products 

(except furniture) - manufacturing 

25. Furniture and fixtures - 
manufacturing 

26. Paper and allied products - 
manufacturing 

27. Printing, publishing, and allied 
industries 

28. Chemicals and allied products - 

manufacturing F 

29. Petroleum refining and related 
industries 


3. 


New York State Inventory 


A. Agriculture 

Areas : 

Ao Orchards 
Av Vineyards 

Ah Horticulture, floricul- 
ture 

Ay Specialty farms 
At High- intensity cropland 
Ac Cropland and cropland 
pasture 

Ap Permanent pasture 
Ai Inactive agricultural 
land 

Ui Other inactive lands 
Uc Lands under construction 


31. Rubber and miscellaneous plastic 
products - manufacturing 

32. Stone, clay and glass products - 
manufacturing 

33. Primary metal industries 

34. Fabricated metal products - 
manufacturing 

35. Professional, scientific, and 
controlling instruments; photo- 
graphic and optical goods; watches 
and clocks - manufacturing 

39. Miscellaneous manufacturing - NEC 


Point data ; 

Ay Specialty farms 
y-1 Mink 

y-2 Pheasant and game 
y-5 Aquatic agriculture 
y-6 Horse farms 
d Dairy farms; number 
e Poultry operation; 
number 

f Active farmsteads; 
number 

Forest Land 

Areas: 


. . Fc Forest brushland 

Fn Forest land 
Fp Plantations 

W. Water Resources 


Areas ; 


Wn Natural ponds and lakes 
(1 acre +) 

Wc Artificial ponds and 
reservoirs (1 acre +) 

Ws Streams and rivers 
(100' +) 

Wh Hudson River 
Km Marine lakes, rivers 
and seas 
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Table A (cont.) 


4. Transportation, communication, 

and utilities 

41. Railroad, rapid rail transit, 
and street railway transpor- 
tation 

42. Motor vehicle transportation 

43. Aircraft transportation 

44. Marine craft transportation 

45. Highway and street right-of- 
way 

46. Automobile parking 

47. Communication 

48. Utilities 

49. Other transportation, communi- 
cation, and utilities, NEC 

5. Trade 

51. Wholesale trade 

52. Retail trade - building mate- 
rials, hardware, and farm 
equipment 

53. Retail trade - general mer- 
chandise 

54. Retail trade - food 

55. Retail trade - automotive, 
marine craft, aircraft, and 
accessories 

56. Retail trade - apparel and 
accessories 

57. Retail trade - furniture, home 
furnishings, and equipment 

58. Retail trade - eating and 
drinking 

59. Other retail trade, NEC 

6. Services 

61. Finance, insurance, and real 
estate 

62. Personal services 

63. Business services 

64. Repair services 

65.. Professional services 

66. Contract construction services 

67. Governmental services 

68. Educational services 

69. Miscellaneous 


Wb Shrub wetlands, bogs, 
marshes 

Ww Wooded wetlands 
Point Data : 

n Natural ponds and lakes; 
number 

c Artificial ponds and 
reservoirs: number 
p Ponds less than 1 acre 
in size: number 
1 Lake shoreline: miles 
s Streams and rivers: 
miles 

N. Nonproductive land 

Ns Sand (unstabilized) 

Nr Rock (exposed) 

R. Residential Land Use 
Areas : 

Rh High density (50' 
frontage) 

Rm Medium density (50-100' 
frontage) 

Rl Low density (100' + 
frontage) 

Re Residential estates 
(5 acres +) 

Rs Strip development . 

Rr Rural hamlet 
Rc Farm labor camp 
Rk Shoreline cottage 
development 

Point data : 

k Shoreline developed in 
cottages: miles 
z High-rise apartment 
buildings: number 
v Trailer parks: number 
x Rural non- farm resi- 
dences never a farm 
residence: number 
o Rural non-farm resi- 
dences once a farm 
residence: number 


85 


Table A (cont.) 


7. Cultural, Entertainment, and C. 

Recreational 

71. Cultural activities 

72. Public assembly 

73. Amusements 

74. Recreational activities 

75. Resorts and group camps 

76. Parks I. 

79. Other cultural, entertainment, 

recreation, and NEC 

8. Resource Production and Extrac- 
tion 

E. 

81. Agriculture 

82. Agricultural related acti- 
vities 

83. Forestry activities and re- 
lated services 

84. Fishing activities and re- 
lated services 

85. Mining activities and related 
services 

89. Other resource production and 
extraction, NEC 

9. Undeveloped Land and Water Areas 

91. Undeveloped and unused land 

area (excluding non- commercial 
forest development) OR. 

92. Noncommercial forest develop- 
ment 

93. Water areas 

94. Vacant floor area 

95. Under construction 

99. Other undeveloped land and 
water areas, NEC 


1 / NEC - Not elsewhere classified 


Commercial Areas 
Areas : 

Cu Central business district 
Cc Shopping center 
Cs Strip development 
Cr Resorts 

Industrial Areas 

Areas : 

Ii Light manufacturing 
Ih Heavy manufacturing 

Extractive Industry 

Areas ; 

Es Stone quarries 
Eg Sand and gravel pits 
Em Metalic mineral extrac- 
tion 

Eu Underground mining 
Point data : 

Eu Underground mining: 
types present 
u-1 Oil and gas 
u-2 Salt 
u-3 Other 
u-4 Abandoned 

Outdoor Recreation 

Areas : 

OR All outdoor recreation 
activities 

Point data : 

OR Outdoor recreation facil- 
ities: types present 
OR-1 Golf courses 
OR-2 Ski areas, other 
winter sports 
OR- 3 Beaches and pools 
OR-4 Marinas, boat launch- 
ing sites 
OR- 5 Campgrounds 
OR- 6 Drive-in theaters, 
race tracks, amuse- 
ment parks 
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Table A (cont.) 


OR- 8 Fairgrounds 
OR-9 Public parks 
OR- 13 Shooting, archery 
OR- 16 Private company 

facilities, community 
areas 

P. Public and Semi-Public Land 
Uses 

Areas : 

P All public & semi-public 
areas 

Point data : 

P Public & Semi-public 
areas - types present 
P-1 Educational institu- 
tions 

P-2 Religious institutions 
P-3 Health institutions 
P-4 Military bases and 
armories 

P-5 Solid waste disposal 
P-6 Cemeteries 
P-7 Water Supply treatment 
P-8 Sewage treatment plants 
P-9 Flood control struc- 
tures 

P-11 Correctional institu- 
tions 

P-1 Road equipment centers 
P-16 Welfare centers, 
county farms 

T. Transportation 

Areas : 

Th Highway interchanges, 
limited access right- 
of-way, etc. 

Tr Railway facilities 
Ta Airport facilities 
Tp Marine port and dock 
facilities 
Ts Shipyards 
Tl Marine locks 
Tt Communication and 
utility facilities 


8? 


Table A (concluded) 


Point data : 

h Highway category: high- 
est present 
h-o None 

h-3 Unimproved, gravel, 
town roads 

h-4 Two- three land highway 
h-5 Four- line highway 
h-6 Divided highway 
h-7 Limited access highway 
h-8 Limited access inter- 
change 

Tr Railway facilities: type 
present 

r-1 Abandoned right-of-way 
r-2 Active track 
r-3 Switching yards 
r-4 Stations and structures 
r-5 Spur 

Ta Airport facilities: type 
present 

a-1 Personal 

a-2 Non- commercial 

a-3 Commercial 

a-4 Airline 

a-5 Military 

a-6 Heliport 

a- 7 Seaplane base 

Tb Barge canal facilities: 
types present 
b-1 Channel 
b-2 Lock 

b-3 Abandoned channel 

Tt Communications & utilities: 
types present 

t-1 TV- radio tower 
t-2 Microwave station 
t-3 Gas & oil-long distance 
transmission 
t-4 Electric power- long 
distance transmission 
t-5 Water - long-distance 
transmission 

t-6 Telephone - long-dis- 
tance transmission 


Table Bl. ouramary of the May test field results showing the number of data points and percent of total 
•data points in each computer class. 
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Total data points for computer classes 











































Summary of the July test field results showing the number of data points and percent of 
total data points in each computer class. 



Total data points for computer classes 































Spectrally separable computer classes 
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able B3. Summary of the September test field results showing the number of data points and percent. of 
total data points in each computer class. 
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appendix c 


Figure Cl. Identification of agricultural cropland 

and an overlay of the land use categories 
present in Flightline 219 . 


Cropland* 


Land use categories** 


Soybeans. . . . . s 
Pasture ..... P 
Idle cropland . . i 

Rye r 

Wheat ...... w 

Hay ....... h 

Corn c 


Row crop ........ A 

Small grain, ...... B 

Hay ...«•••, «*C 
Pasture ........ D 

Forest, continuous . , . E 
Forest, discontinuous . F 

Extraction , G 

Fresh water. ...... H 

Waste water I 

Marshland. ....... J 

Urban core K 

Suburban ........ L 


* Annotated on aerial photograph (Note: penciled numbers and 
letters pertain to field and tract information) 


** Annotated on overlay 
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Figure Cl. 
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Figure Cl . Cont 
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Figure Cl . Concluded 
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